% Run some tests
clear all
close all
clc
% Mallas para los tests
m = 5;
NChen  = [5 9 17 33 65];

% Curvas de convergencia

nu_v = [0.3 , 0.49 , 0.499 , 0.4999 , 0.499999];

for i=1: size (nu_v,2)
    nu = nu_v(i);
    exercise32_FEM;
    L2_conv(i,:) = L2_err;
end
%%
%set(gcf,'Color',[1,0.4,0.6])
legends = char('\nu=0.3','\nu=0.49', '\nu=0.499','\nu=0.4999','\nu=0.499999');
formats = char('-go','-bs','-rv','--mx','--yv','--[0.5 0.5 0.5]v');%, '-.ko', '-.ks');
plotConvergence ( h_n , L2_conv' , legends , formats );